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ABSTRACT 

' The gravitational microlensing light curves that reveal the presence of extrasolar 

■ planets generally yield the planet-star mass ratio and separation in units of the Einstein 

ring radius. The microlensing method does not require the detection of light from the 



00 



planetary host star. This allows the detection of planets orbiting very faint stars, 
but it also makes it difficult to convert the planet-star mass ratio to a value for the 
planet mass. We show that in many cases, the lens stars are readily detectable with 
high resolution space-based follow-up observations in a single passband. When the lens 
star is detected, the lens-source relative proper motion can also be measured, and this 
QQ ' allows the masses of the planet and its host star to be determined and the star-planet 

■ separation can be converted to physical units. Observations in multiple passbands 
provide redundant information, which can be used to confirm this interpretation. For 

. the recently detected super-Earth planet, OGLE-2005-BLG-169Lb, we show that the 

! lens star will definitely be detectable with observations by the Hubble Space Telescope 

(HST) unless it is a stellar remnant. Finally, we show that most planets detected by a 
space-based microlensing survey are likely to orbit host stars that will be detected and 
O ' characterized by the same survey. 

w 

■ Subject headings: gravitational lensing, planetary systems 
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^ . 1. Introduction 

One of the primary strengths of the gravitational microlensing extrasolar planet detection 
method is its sensitivity to low-mass planets at separations of a few AU. This has recently been 
demonstrated with the discovery of two of the lowest mass extrasolar planets known: OGLE-2005 - 



BLG-390Lb, a 5.5^2 7-^^© planet orbiting an M-dwarf in the Galactic bulge ( Beaulieu et al.ll2006l ) 
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and OGLE-20 05-BLG-169Lb, a 13 "^g Mm which probably orbits a K or early M star in the inner 



Galactic disk ([Gould et al.l 120061 ). Both of these planets have a separation from their host star of 
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about ~ 3AU, and this places them well outside the current sensitivity range o f the radial velocity 
method, which is only sensitive to such planets with a separation of ^ 0.7 AU (jRivera et al.ll2005l : 
Lovis et alJbood l. 



The microlensing method is most sensitive to planets at a separation of 1-5 AU, and this 
separation region is particularly interesting from a theoretical point of view because it contains 
the so-called "snow-line" which is an important feature for the core accretion theory of planet 
formation. This "snow-line" is the region of the proto-planetary disk where it is cold enough for 
water-ice to conde nse, and the cor e accretion theory predicts that this is where the most massiv e 
planets will form (jlda Sz LinI 12004 : iLaughlin. Bodenheimer &: Adamd |2004| : iKennedv et al.l |2006| ) . 
According to this theory, giant planets form just outside the "snow line" where they can accrete 
~ IOMq of rock and ice to form a core that grows into a gas giant like Jupiter or Saturn via the run- 
away accretion of Hydrogen and Helium onto this core. However, this theory also predicts that the 
Hydrogen and Helium gas can easily be removed from the proto-planetary disk during the millions 
of years that it takes to build the rock-ice core of a gas-giant. Thus, if the core accretion theory is 
correct, rock- ice planets of ~ lOM^ that failed to grow into gas giants should be quite common . 



Recent core -accretion theory calculations (jLaughlin. Bodenheimer &: Adams 
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Ida Sz Lin 



2004 



Bosal2006al ) predict that it is especially difficult to form gas-giant planets around low-mass stars. 



As of November, 2006, there have been four extrasolar planets discovered by the m icrolensing 
method: two planets more massive than Jupiter (jBond et al.ll2004l : lUdalski et al.ll2005l ) in addition 
to the two ~ lOM^ super-Earth planets. But, the microlensing planet detection efficiency is 
more than an order of magnitude larger for Jupiter-mass planets than for lOM^ planets, so the 
microlensing results i ndicate that planets of lOM^ are significantly more common than Jupiters 



Beaulieu et al 



Gould et al 



200e 



20061 ) and that 16-69% (90% c.l.) of non-binary stars have a planet of ~ lOM^ 
at a separation of ~ 1.5-4 AU. The discovery of these two planets suggests that 
super-Earth planets of ~ lOM,^ are more common than Jupiter mass planets at separations of a 
few AU around the most common stars in our Galaxy. This would seem to confirm a key prediction 
of the core accretion model for planet formation: that Jupi ter mass planets are much more likely to 
form in orbit ar ound G and K stars than around M stars (jLaughlin. Bodenheimer &: Adamsll2004l : 
Ida &: LinI l2004ll . a lthough s uper-Earth planets are also predicted i n the gravitational instability 
model ( Boss 2006b). In fact, ( Laughlin. Bodenheimer &: Adainsll2004) have argued that the Jupiter- 
mass planets found in microlensing events (jBond et al.l l2004l : lUdalski et al.l l2005l ) are more likely 
to orbit white dwarfs than M dwarfs. Clearly, the microlensing detections would provide tighter 
constraints on the theories if the properties of the host stars were known. 

In this paper, we show that the planetary host stars can often be identified in high resolution 
images taken a few years after the microlensing event. In cases where the lens star has been 
identified, we show in §[2]that a complete solution to the microlensing event can generally be found. 
The complete solution gives the lens star and planet masses and the separation in physical units 
instead of just the star-planet mass ratio and separation in Einstein ring radius units. This has 
recently been demonstrated with the first planet discovered by microlensing, OGLE-2003-BLG- 
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235Lb/MOA-2003-BLG-53Lb (|Bennett et alJ 12009 ) . and in § [3l we consider the more favorable 
case of the low-mass planetary microlensing event OGLE-2005-BLG-169, where the lens star will 
be detectable if it is a main sequence sta r of any mass. In § |H we show that a space-based 
microlen sing survev (iBennett &: Rhie II2002I ). like the proposed Microlensing Planet Finder (MPF) 
mission (jBennett et al.ll2004l ) will detect the planetary host stars for most of the detected planets, 
and so, planets detected by MPF will usually come with a complete solution of the lens system. 
We conclude, in § [5l with a discussion of follow-up observations of planetary microlensing events 
and argue that the host stars will generally be detectable unless they are stellar remnants or brown 
dwarfs. 



2. Lens-Source Relative Proper Motion and Complete Solution of Microlensing 

Events 

Gravitational microlensing can be used to study objects, such as brown dwarfs, stellar rem- 
nants, or extrasolar planets that emit very little detectable radiation, but most microlensing events 
provide only a single parameter, the Einstein radius crossing time, Ie, that can constrain the lens 
system mass, distance, and velocity. The situation is significantly improved when the lens-source 
relative proper motion, //reb can be determined, because this yields the angular Einstein radius: 
Ge = ^rei^B- The angular Einstein radius is related to the lens system mass by 

where Di and Ds are the lens and source distances, respectively. Ml = + Mp is the total lens 
system mass: the sum of the star and planet masses. Since -D5 is known (at least approximately), 
eq. [T] can be considered to be a mass-distance relation for the lens star. 

Another constraint on the lens system is needed in order to convert this mass-distance relation 
to a complete solution of the lens event. For most events, this can be accomplished by direct 
detection of the planetary host star. If the brightness of the host (lens) star is measured, then 
the mass-distance relation, eq. [H can be combined with a mass-luminosity relation, such as those 
shown in Figure [H to yield the mass of the lens system. This then allows the individual star and 
planet masses to be determined, since the mass ratio is known from the microlensing light curve. 
This also allows the separation determined from the light curve to be converted from units of the 
Einstein radius. Re, to physical units. 

There is still some residual uncertainty due to the uncertainty in Ds, but this is generally small 
as can be illustrated with some special cases. If <C Ds, then eq. [T] indicates that Ml becomes 
independent of Ds- In the opposite extreme, with the lens very close to the source {Dl ~ Ds), 
we can treat Ds and Ds — D^ as the independent variables and replace D^ hy Ds (except in the 
expression Ds — Di), so in this limit Ml oc D'g/{Ds — Dl) The mass luminosity relations shown 
in Figure [1] can be approximated by a power law over a limited range of masses. So, if we take 
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L oc M" as our mass-luminosity relation, the lens mass is simply given by Ml oc where 
Fi^ is the measured flux of the lens star. (Eq. [T]then gives Ds — Dl in terms of M^.) Since the 
appropriate value of n is almost always n > 2, the fractional error in Ml due to the uncertainty in 
Ds will generally be smaller than the fractional uncertainty in Ds- 

It is important to note that this complete solution of the lens system requires only the lens 
brightness in a single pass-band. Color information obtained by multi-band observations is redun- 
dant, and can be used to constrain the possibility that the source or lens star has a bright binary 
companion. 

An alternative method for the complete solution of a microlensing event is to combine a mea- 
surement of the lens-source proper motion with a measurement of the microl ensing parallax effect. 
This has already been demonstrated f or microlensing events due to single sta rs (jPrake. Cook. &: Keller 
2OO4I 



Gould. Bennett, Alved 12004 ) and stellar binaries (|An et al.ll2002l ). The detection of the 

1 



the 



microlensing parallax effect provides a measurement of fe, the projected Einstei n radius (proj ected 
from the position of the source to that of the observer). As originally shown by iGouldl (|199 
measurement of both 9e and fg yields the lens system mass, 

Ml 



(2) 



without any dependence on Ds or Dl- However, fg i s generally rneasur ed by detecting the ef- 
fect of the Earth's orbital motion in the light curve (jAlcock et al.lll995l ). and this means that 
reliable measurements of fg can generally only be made for very long duration events or events 



where the lens-source motion is unusually small (^Alcqck et al.ll2001bl: iDrake. Cook. &: Keller 
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I2OO5I1. For shorter events, the lieht curve 
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single component of the two-dimensional fe vector. It is then possible to determine the lens mass 
from eq. [2] if the direction of /ij-ei is measured, since these two vectors are parallel. 



2.1. Lens-Source Relative Proper Motion Determination Methods 



There are two methods of measuring the relative proper motion. The first involves a mea- 
surement of finite source effects in the microlensing light curve combined with an estimate of the 
angular size of the source. Th ese finite source effects are commonly detectable in binary lens 
events ( Mao fc Paczynski Ill99ll ). whe n the angular position of the source crosses or comes very 



close to the caustic curve of the lens ( Bennett et al. 



1999 : lAfonso et al.ll2000l : lAlcock et al.ll200d : 



Jaroszvnski et al.ll2005l ). Finite source effects are particularly common in planetary microlensing 
events because the planet is likely to be detected only if the source crosses or closely approaches 



the caustic curve. In fact, all four planetary microlensing event s observed to date ([Bond et al 



2004 



Udalski et al. 



2005 



Beaulieu et al. 



2006 



Gould et al.l I2OO6I ) have revealed finite source ef- 



fects, which reveal the source crossing time, t^,. The angular radius of the source star, 9^, can be 
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estimated from its color and brightness. This yields the relative proper motion, /ifei = 0*/t*, and 
the angular Einstein radius, 9e = O^tE/t^, because tE can virtually always be determined from the 
light curve of planetary microlensing events. 

Even if the light curve does not reveal any finite source effects, it is still possible to measure 
Hrei directly by detecting the lens st ar and measuring th e lens-source separation as was done for 



microlensing event MACHO-LMC-5 (jAlcock et al.ll2001bl ). For this event, the relative proper mo- 



tion was quite large, fi^d = 21.4mas/yr, because the distance to the lens is small, Dl = 550 pc. 
As a result, the lens and source were easily resolved in HST images taken 6.3 years after peak 
magnification. 

Galactic bulge microlensing events with planetary signals generally have relative proper mo- 
tions that are much smaller than this, typically, //rei ~ 5mas/yr. Thus, the lens stars would 
typically not be resolved in HST images taken less than a decade after peak magnification. But, 
due to the stability of the HST point-spread function (PSF), it is possible to measure lens-source 
separations that are much smaller than the width of the PSF. This is accomplished by measuring 
the elongation of the combined lens-source image due to the fact that it consists of two point sources 
instead of one. 



2.2. Image Elongation 



Figure [2] shows a simulated images from the Hubble Space Telescope (HST) Advanced Camera 
for Surveys (ACS) of the source and lens star for microlensing event OGLE-2006-BLG-169 taken 
2.4 years after peak magnification when the lens-source separation is predicted to be 20.2 it 1.7 mas. 
These simulated images represent the co-added dithered exposures from a single orbit of obser- 
vations in the /-band (F814W), and the three different rows of images represent three possible 
masses for the planetary host star. While the 20.2 ± 1.7 mas lens-source separation is smaller than 
the diffraction limited /-band PSF FWHM of ~ 76 mas, the elongation of the PSF of the blended 
source-lens star pair is clearly visible in the residual images shown in the central and right hand 
columns. 

The high precision image elongation measurements that we desire depend on a precise knowl- 
edge of the PSF, but this is to be expected with space-based observations of microlens source 
stars. The Galactic bulge fields where microlensing events are discovered have large numbers of 
relatively bright, but reasonably well separated stars in the field of the lens, and this enables the 
determination of very accurate PSF models, as long as the space telescope observations are dithered 
to correct for image und ersampling. Prior HST programs have dem onstrated that HST has the 



np 



requisite image stability (|Laueijll999l : lAnderson fc King Il200nl . 120041 ) 



There is a potential degeneracy encountered when converting the image elongation into a 
value for the lens-source displacement. The elongated image full-width half-maximum (FWHM) 
in the direction of the lens-source separation is determined by the lens-source separation and their 
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brightness ratio. Thus, if we only measure the increased FWHM of the elongated image, we won't 
be able to determine the brightness ratio and the separation. Fortunately, the microlensing light 
curves with definitive planetary signals will generally provide the information needed to break 
this degeneracy. These light curves have sufficient high precision photometry to determine the 
brightness of the source star from the shape of the microlensing light curve. Furthermore, most 
planetary light curves also exhibit finite source effects, which allow an independent determination of 
the lens-source relative proper motion, /Xi-ei- For these events, we can consider the image elongation 
to be a measure of the the brightness ratio, which can be compared to the source star brightness 
from the light curve as a test for a bright binary companion to the source star. 

The most reliable way to estimate the precision of our measurement of /irei is with Monte Carlo 
simulations using realistic PSF models, which are well understood in the case of HST. However, 
it is instructive two consider a much simpler procedure. Let Ax be the lens-source separation at 
a time, T, after peak magnification. The relative proper motion is therefore /Xrei = Ax/T. If we 
approximate the PSF with a Gaussian profile, g-O-S^^/^Oj then the blended image of the source plus 
lens star will be represented by the sum of two Gaussian profiles with centers separated by Ax. 
We can then approximate this sum of Gaussians with a broader Gaussian, ^-^■^^^h'^ ^ where we 
require that the broader Gaussian has the same RMS value as the original sum of Gaussians (in 
the direction of the lens-source proper motion). This yields 

S^ = sl + fL{l-fL){^xf , (3) 

where fi = Fi/Ftot is the ratio of the flux of the lens to the total flux in the PSF. Since the source 
flux Fs = Ftot — Fl is known from the light curve, and Ftot can be measured from the images being 
analyzed, we can consider to be known from the light curve. The uncertainties in the lens-source 
separation. Ax, and relative proper motion, /ii-ei can be derived. 



Ax /irel V2]VtoI V/l(1-/l)(Ax)2 

Here Ntot is the total number of photons in the PSF. We find that eq. [4] does an excellent job of 
describing the precision of our simulated Ax measurements for a variety of model PSFs, although in 
some cases, it may be necessary to multiply eq. |l]by a fudge factor (to account for the inadequacies 
of our Gaussian approximation). 

This approximate formula is compared to the actual uncertainties, derived from our simulations 
in Figure [3] for parameters appropriate for a s pace-based microlen sing survey, such as the proposed 



Microlensing Planet Finder (MPF) mission (jBennett et al.ll2004l ). (Such a mission would have a 



PSF FWHM 2-3 times worse than HST, but would have the ability to combine images with a 
combined exposure time of months.) 

We expect that eq. [3]may fail to describe the measurement uncertainty for Ax when the lens- 
source separation grows to Ax ^ sq and the lens and source images begin to be separately resolved 
in the high resolution images. In this case, there are more features of the image that can aid in the 
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determination of Ax, so we should be able to measure Ax even better than eq.H] predicts. However, 
this is only likely to occur in cases where ctax is already quite small, and so eq. H] will simply give 
a very conservative estimate of ctax in these cases. 

For most planetary events, the light curve shows finite source effects that allow t^, and therefore 
/^rei and Ax to be determined from the light curve. In these cases, the image elongation can be 
used to measure /l, and the uncertainty in the fi determination is given by 



In most cases, this determination of Jl will be redundant with the determination from the light 
curve fit and the total lens plus source flux. However, if the source star has a bright binary 
companion, at a separation > 10 so that the light curve is not affected, then these two methods 
will yield different results. The binary companion will contribute to the Fs = i^tot(l — fi) value 
as measured by image elongation, but it will not contribute to the value of Fs determined from 
the light curve. Thus, the different measures of Fs provide a measurement of the brightness of the 
source star's binary companion. 

In order to determine how precisely we can measure both Ax and fi from the blended image 
of the source plus lens, we will need to measure the third moment, or skewness, of the blended 
image profile. The third moment is given by 

Q = ^Y.(^-^of , (6) 

where the sum is over the positions of the detected photons and xq is the mean position of the 
detected photons. If the image profile is approximately a Gaussian, then Q can be measured to a 
precision of cjq = s'^ ^/Tb/Ntot- We can also derive the following expression for Q in terms of /l 
and Ax, 

Q = /l(1-/l)(1-2/l)(Ax)3 . (7) 

This expression can be used to derive the following expressions for the uncertainties in and Ax 
when we solve for both, 



(TAx _ / 15 1 / g \3 
Ax "^yNtoth VAx/ 



(8) 



and 



60 / s 



Ntot VAx/ 

These expressions assume that /l ^ 1 and that the uncertainties are dominated by our ability to 
measure Q, which will generally be the case when Ax/s ^ 1. 

An additional complication is the possibility of blending by stars that played no role in the 
microlensing event, which happen by chance to lie close to the positions of the source and lens stars. 
However, the mean separation of stars in the Galactic bulge fields that are searched for microlensing 
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events is in the range of 0.3-1.0 arc seconds, which is much larger than the lens-source separations. 
Thus, if the lens-source star separation can be measured, it is unlikely that stars unrelated to the 
microlensing event will interfere significantly with the lens-source separation measurement. The 
only exception would be stars that are physically associated with the lens or source stars, and this 
possibility is dealt with in ^ l2.5l below. For a space-based microlensing survey, it should generally be 
possible to directly detect the unre lated blended stars w ith a point-source decomposition algorithm. 
(This is similar to deconvolution (iMagain et al.ll2006l ) except that it makes use of the fact these 
Galactic bulge field contain very few sources that are not point-like.) 

We should also note that this estimate of our ability to measu re the separation o f two unresolved 
stars is substantially better than the conservative estimate of iHan et al.l (|2006l ). Our estimate 
i s justified by our detail ed simulations, which are, in turn, based on experience with HST data 



Anderson &: King II2004I ). Thus, we are confident that the conservative assumptions of iHan et al 
200d) are not needed. 



2.3. Color Dependent Centroid Offset 

Another method to detect the effects of the relative lens-source proper motion relies upon the 
likelihood that the lens and source will have different colors. This means that the lens star will 
contribute a different fraction to the total light seen in the source plus lens image blend as seen in 
different passbands. This magnitude of the observed color-dependent image center shift for the V 
and / passbands is 

Axv^i = if I - fv) Ax . (10) 

Axv-i is directly measurable from high resolution HST images, but // and fv depend on bright- 
nesses of the lens and source in both passbands. The brightness of the source star in multiple 
passbands can generally be obtained by multicolor observations at several different magnifications 
during the course of the microlensing event, but the brightness of the lens star is an unknown that 
must be solved for along with Ax. Equations [T] and [10] can be combined with the V and / band 
mass-lumonisity relations shown in Figure [1] to yield four equations for the four unknowns: M, 
Vl, /lens 5 and Ax = OET/tE- These equations can then be solved to yield a complete solution for 
the microlensing event par ameters. This effect has been detected for the first extrasolar planet 



discovered by microlensing (jBennett et al.ll2006l ). but in this case, firei (and therefore 9e and Ax) 



was already known from the microlensing event light curve. 

An important advantage of this method is that for small Ax, the signal is oc Ax. Thus, as 
long as the colors of the source and lens stars differ significantly, this method should measure small 
Ax vlaues more precisely than the image elongation method. 
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2.4. Interstellar Extinction 



Because our methods involve measuring the brightness and color of the planetary host star in 
the direction of the Galactic bulge, interstellar extinction is an important consideration. Extinction 
affects the determination of the brightness and color of both the source and lens stars. The angular 
radius of the source star, 0*, is determined from its brightness and color, and 9^ is used to determine 
the lens-source relative proper motion, /ij-ei, and the angular Einstein radius, 6e- Fortunately, 
extinction makes stars appear both fainter and redder, and these two effects push our estimate of 
in opposite directions. As a result, our estimate of 9^ is not highly dependent on the uncertainty 
in our extinction estimate. 

The dust in the Gal actic disk is generally modeled with a disk scale height of ~ 100 pc 
y). Since the source and lens stars are generally located at a distance 



(jPrimmel fc Spergel 



2001 



^ 1-2 kpc, and are ^ 2° from the Galactic plane, the extinction is primarily in the foreground of 
both the source and lens star. The extinction toward the Galactic bulge is also quite patchy, so it 
cannot be reliably estimated by a simple model. Instead, the standard practice for the interpreta- 
tion of microlensing events is the estimate the extinction toward the source based upon the observed - 
colors and magnitudes of stars within 1- 2 arc minutes of th e microlensed source (jBond et al.ll2004l : 
Udalski et al.ll2005l : iBeauheu et al.ll2006l : iGould et al.ll2006l ). For the lens star, the extinction must 
be modeled with a probability distribution, which mat ches the measured extinct ion at the dis- 
tance of the source and follows the exponential scaling of lDrimmel Spergell (|200lh . The variance 
of the probability di stribution accounts for the patchiness of the actual distribution. In practice 
(jBennett et al.l |2006| ) , this procedure amounts to only a slight modification to the simpler case 
in which all the extinction is assumed to lie in the foreground of both the lens and source. For 
simplicity, in this paper, we will assume that extinction for the lens is identical to the extinction 
for the source star, and we will not explicitly discuss the extinction in the remainder of this paper. 



2.5. Binary Companions to the Lens or Source 



Although binary stars represent a minority of star systems (jLadal |2006| ) , it is important to 
consider the effect of a binary source or lens star on our analysis of the planetary host star properties. 
A binary companion to the lens or source star would generally be unresolved from the source, so it 
would complicate our analysis. But, such a possibility is also constrained by the properties of the 
microlensing light curve. 

A binary companion to the source star can often be detected via light cur ve oscillations due 



to orbital motion o f the source star if the orbital period is less than a year (jPerue et al. 



1999 



Alcock et al.ll2001al ). Microlensing events with main sequence source stars and detected planetary 
light curve features generally have a peak magnification of j4max ~ 10 because lower magnification 
events typically have poor photometry due to the crowding in ground-based images. If the event 
is well sampled near peak magnification, binary source effects in the light curve will be visible 
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with an amplitude of ~ a/(i?£;^niax)) where a is the semi-major axis of the source star orbit. For 
^max > 10, this imphes that even the orbital motion of a contact binary will be visible if the 
photometric precision is close to 1%. Thus, in most cases, a binary companion to the source must 
have a semi-major axis ^ 1 AU in order to avoid visible light curve effects. If the source star's 
companion is not much fainter than the source, then the lensing of the companion star may be 
visible in the light curve unless the source-companion separation is ^ Re ^ several AU. If the 
separation of the source from its companion is ^ 5-10 mas, depending on its brightness, then we 
will be able distinguish its position from that of the source star using the same methods that we 
have discussed to detect the lens-source offset. For a typical bulge source distance of 9kpc, this 
corresponds to a source-companion separation of 45-90 AU, so for separations larger than this, the 
binary companion will affect our measurements of lens-source astrometry. 

The light curve constraints on a possible companion to the lens and planetary host star depend 
more sensitively on the parameters of the specific microlensing event. For a very high magnification 
event, such as OGLE-2005-BLG-169 (see §[3]), a companion of similar mass to the planetary host 
star can be excluded over a very wide range of separations from ~ 0.01-lOOii^; or ~ 0.03-300 AU. 
For events of more modest magnification, the binary separations are excluded over a more modest 
range of separations, ~ 0.1-lOi?^;. So, there will usually be at least a slim chance of the lens star 
system is a close binary star system with a semi-major axis o <C 1 AU or a wide binary system 
with o ^ 1 AU. The observational signature of a wide binary companion is similar to the signature 
of a binary companion to the source in that there is an extra source of light located very close to 
the positions of the lens and source. However, this extra source will move with the planetary host 
star, and so the image elongation signal will be stronger. In contrast, if the planetary host star is 
a close binary, then both stars will contribute the lens mass in eq. [21 so the lens star system will 
be fainter and redder than implied by the single-star mass-luminosity functions shown in Figure [TJ 



2.6. Full Solution with Combined High Resolution Follow-up Data 

One complication that may tend to interfere with the recovery of the planetary host star 
parameters is the possibility that either the source star or the planetary host star will have a 
binary companion star that contributes a significant amount of flux to the blended image. Such a 
star would generally be unresolved from the position of the lens and source star, and this would 
hamper our ability to recover the brightness of the lens star. Generally, such a binary companion 
to the source or lens star would have to have a separation ^ 10 AU from the lens or source to avoid 
an obvious light curve signal. For high magnification microlensing events, a binary companion 
would have to have an even larger separation to avoid detection, but for a low magnification event, 
it would be possible to find a planet orbiting a very close binary star system with no detectable 
light curve features. Fortunately, the high resolution follow-up data provides more than enough 
constraints to identify these binary companions, in most cases. 

The critical unknown parameters for a microlensing event are the lens mass, distance and rela- 
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tive proper motion, fJ^j-eU and these are constrained by the hght curve parameters Ie and (which 
is usuahy, but not always, measured). The hght curve also yields the source star brightness, usually 
in two different pass bands, but these measurements can only constrain the lens star parameters 
in conjunction with the high resolution images. With deep, two-color high resolution space-based 
follow-up observations, it is possible to measure five additional parameters: the source plus lens 
brightness in two colors, the image elongation in two colors, and the source plus lens centroid offset 
between the two colors. In all, we have 8-9 constraints for the 3 microlensing event unknowns, so 
the problem is significantly over-constrained. 

These 5-6 additional parameters can be used to constrain other possibilities, such as a bright 
binary companion to the lens or source star. Such a companion would add 1 parameter (the 
companion mass) to the problem or possibly three parameters if the separation of the companion 
is large. We could probably constrain a companion to the lens, source or both with these 5-6 
additional parameters. However, the best way to solve a lens system that involves a bright, distant 
lens or source companion would be to obtain high resolution space-based images at multiple epochs. 
This would add four additional measured parameters and allow us to distinguish stellar separation 
due to the lens-source relative proper motion from the separation of a distant binary companion to 
the lens or source. 

Of course, we haven't shown that these parameters will always be measured to sufficient 
precision to guarantee a complete lens solution. A complete investigation of this question is beyond 
the scope of the present paper, but we will demonstrate that the question of binaries is easily 
resolved for the example event that we discuss in § [3l Also, because Ax grows with time, the 
precision of our measurements will grow with time, so it is likely that many ambiguities can be 
resolved with additional follow-up observations. Thus, we can generally expect a complete solution 
of the microlensing event when light from the lens can be observed by HST or another high angular 
resolution spaced-based telescope. 

In practice, some of these measurements can have significant or correlated uncertainties, so the 
observational constraints are best applied through a Bayesian analysis. In fact, such an analysis is 
able to provide relatively precise event parameters even without even observational constraints to 
constrain binary companions to the source and lens, because the a priori probability of a binary 
companion of similar brightness t o the source is relativ ely small. For example, in the case of the 



first microlensing planet discovery iBennett et al.l (|2006l ) were able to constrain the host star mass 
to an accuracy of about 13 %. 



3. A Low-Mass Planetary Event w^ith a Detectable Lens Star 

Three of the four planets discovered by microlensing to date had main sequence source stars, 
which enables detection of the lens star with HST images a few years after the event. Of these three 
planets, OGLE-2005-BLG-169Lb is the most interesting, since it is a cool, Super-Earth mass planet 
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with Mp = 13_gM0. Such a planet would be invisible to other planet detection methods, and the 
light curve analysis indicates that it does not have a gas giant companion in the separation range 
1-10 AU. Statistical arguments indicate that such cool, Su per-Earths are the most common type of 
extrasolar planet yet to be discovered (jGould et alj 120061 ) . However the detailed properties of this 
planet and its host star are uncertain because the host star has not been detected. In the absence 
of the host star detection, the probability distributions for the star and planet masses, distance, 
and separation can only be determined by a Bayesian analysis. This analysis makes use of the 
parameters from the microlensing light curve, including the lens-source relative proper motion of 
//rei = 8.4 lb 0.6mas/yr, and the results are show n in Figured These have assumed a Han- Gould 
model for the Galactic bar (jHan &: Gouldlll995l ). a double-exponential disk with a scale height of 
325 pc, and a sca l e leng th of 3.5 kpc, as well as other Galactic model parameters as described in 
Bennett &: Rhie (120021 '). Because this model is slightly different from the Galactic model used 



bv iGould et al.l (|2006l ). the resulting parameters differ slightly from their results. We find a lens 
system distance of Dl = 2.7^5'p kpc, a three dimensional star-planet separation of a = 3.3^i'q AU 



and main sequence stellar and planetary masses of = O.52|^q22M0 and Mj 
assume that white dwarfs have an a priori probability to host planets that is equal to that of main 
sequence stars (at the separations probed by microlensing), then there is a 35% probability that 
the host star is a white dwarf. The possibility o f a brown dwarf h ost star is excluded by the light 
curve limits on the microlensing parallax effect ([Gould et al.l 120061 ). 
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Figure Sl^d) shows the probability distribution of the /-band magnitude of the planetary host 
star compared to the source star at / = 20.58 ib 0.10. The implied planetary host star brightness 
distribution has a median and l-o" range of I\ens = 21.9 , but the most interesting feature of this 
figure is that the probability of a main sequence lens fainter than / = 23 vanishes. This is because 
the mass-distance relation, eq. [1] ensures that the lens star will be nearby and at least at bright as 
/ = 23, even if it is at the bottom of the main sequence at M^, = O.O8M0. In fact, the microlensing 
parallax constraint from the light curve yields a lower limit for the lens star mass of M^, ^ O.14M0. 
Thus, as the upper panel of Figure [5] indicates, the planetary host star must be at least 16% of 
the brightness of the combined lens plus source star blended image, and this implies that it will 
be detectable if it is not a stellar remnant. So, for the faintest possible main sequence lens star 
of the minimum mass, M^, = O.14M0, eq. [5] implies that /l can be determined with a precision 
of 0.017 in a single orbit of HST ACS/HRC observations in the F814W pass band (with 97,000 
detected photons). This allows a tight constraint on a bright binary companion to the source star. 
The situation is a bit worse for ~ 0.5 when the source and lens are almost equally bright, but 
in this case. Figure [5] shows that the lens star will be easily detectable in the V and B passbands, 
where the HST PSF (sq) is considerably sharper. So, we don't anticipate a significant problem in 
measuring f^, as long as we have a measurement of Fs in more than one passband. 

If we assume the most likely case of no detectable companion to the source, then eq. |4] implies 
that //rei can be measured with a precision of 4.3% in a single orbit of HST ACS/HRC observations 
in the F814W pass band (with 97,000 detected photons). This is an improvement over the 7% 
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measurement of fi^ei that comes from the microlensing hght curve, and would provide independent 
confirmation of the planetary interpretation of the light curve. 

The color dependent effects visible in high resolution images of the OGLE-2005-BLG-169 lens 
and source star blend are summarized in Figure [Sj which is based upon the measured relative lens- 
source proper motion for this event. The top panel shows the fraction fi of the source plus lens 
light that is contributed by the lens for different HST ACS/HRC passbands, and the bottom panel 
shows the predicted color dependent image center shifts. Both of these are shown as a function of 
the planetary host star mass, and the dashed grey lines indicate the constraints on the host star 
mass, 0.14Mq < < 0.76Mq, implied by the existing observational limits on the microlensing 
parallax effect and the maximum brightness of the lens star. 

Several features are apparent in Figure El First, a measurement of a color dependent centroid 
shift does not yield a unique mass, because there are usually two different masses that yield the 
same centroid shift. This occurs because the centroid shift can become small if the lens star is faint 
or if the lens star has a similar color to the source star. However, this will generally not create 
an ambiguity in the interpretation of events because the brightness of the source star is known. 
The "faint lens" solution implies a fainter total brightness for the lens plus source blend than the 
"similar color" solution. Thus, the centroid shift degeneracy is broken by the constraint on the lens 
star brightness, which we get from the light curve. Also, if the system is observed in more than 
two passbands, the multiple color dependent centroid shifts will also resolve this degeneracy. 

Figure [5] also indicates that it will be difficult to determine the OGLE-2005-BLG-169L lens 
star mass precisely if it is in the range 0.2Mq ^ M^, ^ OAMq, because in this mass range the 
increase in the /-band lens brightness with distance due to the mass-distance relation (eq. [1]) is 
nearly compensated for by the decrease in brightness due to the greater lens distance. This effect 
is also responsible for the peak at I\ens — 21.9 in Figure W^d). This ambiguity can be resolved 
with observations in other passbands. The V and i?-band lens brightness fractions and the B — V 
centroid shift curve does not have this ambiguity. However, the lens is also predicted to be rather 
faint in these bands if the lens mass is in the ambiguous range. So, it would be sensible to attack this 
question with a later set of observations if the lens star appears to lie in the Q.2Mq ^ ^ OAMq 
mass range to take advantage of the larger lens-source separation at later times. Note that this 
near degeneracy does not occur for planetary host stars that reside in the Galactic bulge because 
the mass-distance relation becomes quite steep when approaches Ds- 

4. Complete Lens Solutions from a Space-Based Microlensing Survey 

While follow-up observations with HST or JWST may be feasible for a handful of planets 
discovered each year by microlensing, it will probably be difficult to obtain enough observing time 
to follow-up all the planets discovered by microlensing when the discovery rate increases. Certainly, 
with the expected discovery rate of many hundreds or even a thousand planets per year from a 



- 14 - 



dedicate d space mission (IBennett fc Rhie 1 12002| ) . like the proposed Microlensing Planet Finder 



(MPF) ([Bennett et al.ll2004l ). it seems virtually certain that only a very small fraction of events 
could be followed up with observations with HST or JWST. Fortunately, a space-based microlensing 
survey will have the capability to perform its own follow-up observations. 

The angular resolution of a dedicated space-based microlensing survey, like MPF, will be worse 
than HST's resolution due to its smaller aperture (1.1m vs. 2.4m) and near-IR passband. However, 
the decreased angular resolution is more than compensated for by the much greater total exposure 
time that will be provided by a dedicated space telescope. While the MPF PSF may be three times 
the width of the HST PSF, but the number of detected photons in a 1-month MPF exposure can 
be a factor of 3000 times the number of photons detected in a single HST orbit. For fixed fi, eq. H] 
gives 

^oc^JL— (11) 

SO we expect that MPF will be able to measure /irei about 6 times more precisely than a single 
orbit HST observation (for fixed Ax). 

A space-based survey, like MPF, will have several methods that can be used to obtain the 
complete solution of the microlensing event. These include: 

1. Combining the mass-distance relation, eq. [1] with a main sequence mass-luminosity relation 
similar to those shown in Figure [H as described in Section §[221 usually gives the most precise 
results. The angular Einstein radius, Be can be determined from the measurement of finite 
source effects in the light curve {i.e. by determining t*), or by a measurement of the image 
elongation, as described in 



The masses and distances can also be estimated from the brightness and color measurements 
of the lens light as seen by the space-based survey. This method is not generally as precise 
as method 1, but it can be considered to be an independent check on method 1. It can also 
constrain the possibility of a binary companion to the lens star that is too close or distant 



np 

for its lensing effects to be apparent ([Han[ [20051 ). as discussed in ^ 12.51 and 



3. The lens masses can be determined by combining measurements of /Urei and the microlensing 
parallax effect. This has the advantage that it can be used even when the lens star is too 
faint to detect, but a complete microlensing parallax measurement will only be available for 
events that are longer than average. 

The two main methods for determining the lens-source relative proper motion, ^reh differ in 
that the measurement of the source radius crossing time, i*, only determines the magnitude of 
/ireh while the measurement of the image elongation determines the full two dimensional vector 



/Xroi- This is useful because it is generally much easier to measure only a single component ([Gould 



19981 ) of the projected Einstein radius, fe, which is a two-dimensional vector with a length equal to 



the Einstein radius projected to the position of the Sun and a direction parallel to the transverse 
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component of the lens-source relative velocity. However, since \\ /^reb the measurement of a 
single component of fe can be combined with the direction of /^treb yield the magnitude of the 
projected Einstein radius, which is needed to determine the lens system masses via eq. [2l 



We have simulated (iBennett &: Rhie I l2002l ) the planetary events expected for a space-based 



microlensing survey similar to the MPF mission (jBennett et al.ll2004l ) to determine how well the 
parameter of the detected planetary systems can be determined. The results of these simulations 
are displayed in Figure [6l The parameters are solved for using methods 1 and 3 above, and the 
method that gives the most precise parameters is reported in Figure [6j One simplifying assumption 
that we have made is to assume that only the single component microlensing parallax asymmetry 
can be measured, instead of the full fg vector. This means that we don't have a direct method to 
determine the lens masses for planetary host stars that are too faint to detect. (These are generally 
white dwarfs, brown dwarfs, and late M-dwarfs). If we assume that planetary masses scale in 
proportion to the host star mass, but does not otherwise depend on stelar type, then about 24% 
of stars would fall into this "invisible star" category. 

Figure [6] considers only planets found orbiting main sequence host stars, and panel (a) in- 
dicates that the distribution of host star masses is rather flat with precise mass determinations 
most common for the more massive host stars. Panels (b) and (c) show the distributions of the 
separation and mass uncertainties from these simulations with median uncertainties of 5.2% and 
10.2%, respectively. These simulations clearly show that for the majority of planets discovered by 
a space-based microlensing survey, the star and planet masses, separation and host star type will 
be determined with reasonable precision. 



Discussion and Conclusions 



We have shown that the main uncertainties in the properties of planetary microlensing events 
can be overcome with space-based observations. For planets detected with ground-based microlens- 
ing observations, space-based follow-up imaging can detect the lens star, which allows the properties 
of the star and its planet to be determined. For a space-based microlensing survey, no follow-up 
observations are needed because the survey data will contain the information needed to identify 
the lens star and solve for the masses and separation. 

The space-based observations in this paper apply primarily to events with main sequence 
source stars. Main se quence stars are the prime targets for space-based micr olensing surveys 
(IBennett &: Rhie II2002I') and for ground-based attempts to find earth-mass planets (jBennett &: Rhie 



1996 



Wambsgansslll997l ). However, giant source stars do allow the detection of planets a bit more 



massive than the Earth, inc luding the lowest mass planet detected to date, OGLE-2005-BLG- 
390Lb, (IBeaulieu et al.l 120061 ). around a main sequence star. (Of course, much lower rnass p lanets 
have been discovered orbiting a pulsar (IWolszczan &: Fraillll992l : iKonacki &: WolszczanI 120031 ).) 



While giant source stars are attractive targets for microlensing planet searches because their 
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photometry is less affected by blending and because they require much shorter exposure times for 
precise photometry, they have the drawback that the source star is generally brighter than the lens 
star by a large factor. In eq. HI if ^ F/, then Ntot oc Fg and /l oc l/Fs-, so (T^//irei l/\/T^. 
Thus, since typical red clump giant sources are brighter than the typical main sequence source 
stars by a factor of ~ 100, we expect that firei could be measured about ten times worse for a 
giant source than for the same event parameters with a main sequence source star. It would 100 
times more exposure time to compensate for this, which seems implausible in all but the most 
favorable cases. A more reasonable strategy would be to wait longer to allow Ax to grow, or to 
switch to an instrument with muc h higher angular res olution than HST, such as the Very Large 
Telescope Interferometer (VLTI) (jScholler et al.l l2006l). While curr ent VLTI instrumentation is 
unable to observe the lens-source separation (jPelplancke et al.ll200ll ). it is expected that a future 



VLTI instrument will be capable of detecting the OGLE-2005-BLG-390L lens star sometime in 
the next decade (Beaulieu et al., in preparation). It is also possible that high resolution adaptive 
optics imaging could be used to make some of the measurements that we have described here for 
events with main sequence source stars. This would require either very precise knowledge of the AO 
imaging PSF, or a longer delay between the planetary event and the follow-up imaging to provide 
a larger lens-source separation. 

Our investigation has been incomplete in that we have only considered a subset of the infor- 
mation that might come from analysis of the microlensing parallax effect. For a small subset of 
very long timescale events, it is likely that the lens system parameters can be characterized much 
more precisely than we have estimated here. It is also likely that many events with a parallax 
asymmetry measurement, but without a measurement of the two-dimensional fe vector will be 
able to have their parameters measured determined due to future high precision measurements of 
/ij,gi II fe w ith future high angu lar resolution instruments like the VLTI or the James Webb Space 
Telescope (iGardner et al.l 120061 ) with the added benefit of a long time baseline. 



This paper has focused on the situation in which signals of both the lens star and planet are 
detected, but it is also possible to detect planets tha t are so far from their pa r ent star that th e 
star does not yield a photometric microlensing signal (IDi Stefano &: Scalzdll999l : iHan et al.ll2005l ). 
In such a situation, it may be difficult to distinguish the light curves of planets in very distant 
orbits from those of unbound planets that have been ejected from the planetary systems that they 
were born in. But in the former case, the planetary host star will often be detectable by the 
methods described in this paper. If these methods should fail to detect the host star, follow-up 
observations with the Space Interferometry Mission (SIM) will allow the host star to be discovered 
via its astrometric gravitational microlensing sig nal (|ianll200(j l. which will work even for dark host 
stars, such as white or brown dwarfs. However, SIM observations of main sequence source stars 
may be limited by the long exposure times required for sources fainter than V = 19. 

It is often said that microlensing planet discoveries cannot be followed up, but this is not 
completely true. While a repeat of the photometric planetary signal will not usually occur for 
another ~ 10^ years, we have shown that most planetary events provide a prediction of the lens- 
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source relative proper motion, ^Ui-ei, which can be confirmed with high angular resolution follow-up 
observations with HST. Furthermore, these follow-up observations also allow the complete solution 
of the microlensing event, which includes the conversion of the planetary mass ratio into physical 
masses for the planet and its host star, as well as the projected star-planet separation in physical 
units. Finally, we have also shown that a dedicated space-based microlensing survey, such as MPF, 
will collect the data necessary to extract these follow-up observations in its own photometric survey 
images. Thus, a space-based microlensing survey will provide full lens event solutions, with planet 
and star masses, their separation in physical units, and their distance from us, for most of the 
planets that are discovered. 
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Fig. 1. — The adopted mass-luminosity relations in the I (F814W), V (F555W ), and B (F435W) 
passbands. These an empirical mass-luminosity rela,tions (iKroupa &: Toutlll997l) exte nded to highe r 
masses and the i?-band using iDrilling fc Landoltl (|200d ) , ISchmidt-Kalen 1 19821 ) and iGravl (|l992l ) . 
Approximately Solar metallically is appropriate for the lens stars in the inner Galactic disk and 
the bulge. 
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Fig. 2. — Simulated image stacks of multiple dithered exposures of the OGLE-2005-BLG-169 source 
and lens star 2.4 years after peak magnification using the HST/ACS High Resolution Camera (HRC) 
in the F814W filter band. The top row of images assumed a host star mass of M* = O.OSMq, the 
middle row assumes = O.SSMq, and the bottom row assumes M^, = O.GSM©. In each row, the 
image on the left shows the raw image stack sampled at one half the native HRC (28 mas) pixel 
size. The central column shows the residuals after subtraction of the best fit PSF model, showing 
the blended image elongation along the x-axis due to the lens-source separation. The right hand 
column shows these residuals rebinned to the 28 mas native pixel scale. 
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Fig. 3. — The uncertainty in the lens-source separation, Ax, is shown as a function of Ax. The 
different colors indicate different fractions, /l of the blended lens plus source light that is due to 
the lens. Solid lines indicate the results of our simulations, and dashed lines are calculated from 
eq. m These estimates assume a total of 10*^ detected photons from the blended lens plus source. 
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Fig. 4. — Bayesian probability densities for the properties of the planet and its host star if it is a 
main sequence star, (a) The masses of the lens star and its planet (M* and Mp respectively), (b) 

the separation, (c) their distance from the observer {Di); and (d) the I-band brightness of the host 
star. The dashed vertical lines indicate the medians, and the shading indicates the central 68.3% 
and 95.4% confidence intervals. All estimates follow from a Bayesian analysis assuming a standard 
model for the disk and bulge population of the Milky Way, the stellar mass function of Bennett &: 
Rhie (2002). 
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Fig. 5. — The top panel shows the fraction of the source+lens flux that is predicted to come from 
the lens in the I (F814W), V (F555W), and B (F435W) passbands as a function of lens mass. The 
bottom panel shows the predicted color-dependent centroid shifts as a function of mass for 2.4 years 
of relative proper motion at fi^ei = 8.4mas/yr. The grey dashed lines indicate the upper and lower 
limits on the lens mass due to the upper limits on the microlensing parallax and lens brightness. 
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Fig. 6. — (a) The simulated distribution of stellar masses for stars with detected terrestrial planets. 
The red histogram indicates the subset of this distribution for which the masses can be determined 
to better than 20%. (b) The distribution of uncertainties in the projected star-planet separation, 
(c) The distribution of uncertainties in the star and planet masses. Note that it is the two- 
dimensional projected separation that is presented here, and we have not included the uncertainty 
in the separation along the line-of-sight as was done in Figure HI 



